Isospin dependent properties of asymmetric 

nuclear matter 



P. Roy Chowdhury^'\ D.N. Basu^'^ and C. Samanta ^''^''^'^ 

^Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata 700 064, India 

^Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700 064, India 

"-^ Physics Department, Virginia Commonwealth University, Richmond, VA 

23284-2000, U.S.A. 

Physics Department, University of Richmond, VA 23173, U.S.A. 
Abstract 

The density dependence of nuclear symmetry energy is determined from a system- 
atic study of the isospin dependent bulk properties of asymmetric nuclear matter 
using the isoscalar and the isovector components of density dependent M3Y in- 
teraction. The incompressibility for the symmetric nuclear matter, the isospin 
dependent part Kasy of the isobaric incompressibility and the slope L are all in 
excellent agreement with the constraints recently extracted from measured isotopic 
dependence of the giant monopole resonances in even-A Sn isotopes, from the neu- 
tron skin thickness of nuclei and from analyses of experimental data on isospin 
diffusion and isotopic scaling in intermediate energy heavy-ion collisions. This work 
provides a fundamental basis for the understanding of nuclear matter under extreme 
conditions, and validates the important empirical constraints obtained from recent 
experimental data. 
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1 Introduction 



Measurements of nuclear masses, densities and collective excitations have al- 
lowed to resolve some of the basic features of the equation of state (EoS) of 
nuclear matter. However, the symmetry properties of the EoS due to differ- 
ing neutron and proton numbers remain more elusive to date and study of the 
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isospin dependent properties of asymmetric nuclear matter and the density de- 
pendence of the nuclear symmetry energy (NSE) remain the prime objective 
[1,2,3,4,5,6,7,8,9,10]. Many radioactive ion beam facilities and their upgrades 
already exist and many more are being constructed or planned, including the 
Facility for Rare Isotope Beams (FRIB) in U.S.A. [11], SPIRAL2/GANIL in 
France [12], FAIR/GSI in Germany [13], Radioactive Ion Beam (RIB) Fac- 
tory at RIKEN in Japan [14] and the Cooling Storage Ring (GSR) facility at 
HIRFL in Ghina [15]. These new facilities provide the possibility of explor- 
ing the properties of nuclear matter and nuclei under the extreme condition 
of large isospin asymmetry. Gonsequently, the study of the isospin degree of 
freedom in nuclear physics has recently received a lot of importance and atten- 
tion. The ultimate goal of such study is to extract information on the isospin 
dependence of in-medium nuclear effective interactions as well as the EoS of 
isospin asymmetric nuclear matter, particularly its isospin-dependent term or 
the density dependence of the nuclear symmetry energy. This knowledge, es- 
pecially the latter, is important for understanding not only the structure of 
radioactive nuclei, the reaction dynamics induced by rare isotopes, and the 
liquid-gas phase transition in asymmetric nuclear matter, but also many criti- 
cal issues in astrophysics [3,4,16,17]. With recent developments in constraining 
the isospin dependent properties of asymmetric nuclear matter, especially the 
density dependence of the nuclear symmetry energy, it is of great interest to 
see to what extent the results from a theoretical model are consistent with 
these constraints. 

In this work, based on the theoretical description of nuclear matter using the 

density dependent M3Y-Reid-Elhott effective interaction [18,19] (DDM3Y), 
we carry out a systematic study of the isospin-dependent bulk properties of 
asymmetric nuclear matter. In particular, we study the nuclear incompressibil- 
ity Kfx, for symmetric nuclear matter (SNM), the density dependence of the 
nuclear symmetry energy and extract the slope L and the curvature Ksym pa- 
rameters of the nuclear symmetry energy and the isospin dependent part Kasy 
of the isobaric incompressibility. We compare the results with the constraints 
recently extracted from analyses of the isospin diffusion data from heavy-ion 
collisions based on the isospin and momentum-dependent IBUU04 transport 
model with in-medium nucleon-nucleon (NN) cross sections [9,20,21], isoscal- 
ing analyses of isotope ratios in intermediate energy heavy-ion collisions [22] 
and measured isotopic dependence of the giant monopole resonances (GMR) 
in even-A Sn isotopes [23] and from the neutron skin thickness of nuclei [24] . 



2 The nucleeir equation of state and symmetry energy 

The nuclear matter EoS is calculated using the isoscalar and the isovector 
components [25] of M3Y interaction along with density dependence. The den- 
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sity dependence of the effective interaction, DDM3Y, is completely determined 
from nuclear matter calculations. The equilibrium density of the nuclear mat- 
ter is determined by minimizing the energy per nucleon. The energy variation 
of the zero range potential is treated accurately by allowing it to vary freely 
with the kinetic energy part e'^*" of the energy per nucleon e over the entire 
range of e [26]. This is not only more plausible, but also yields excellent re- 
sult for the incompressibility K^o of the SNM which does not suffer from the 
super luminosity problem. 

Assuming interacting Fermi gas of neutrons and protons, with isospin asym- 
metry X = y^^, p — Pn + Pp, where p„, pp and p are the neutron, 
proton and nucleonic densities respectively, the energy per nucleon for isospin 
asymmetric nuclear matter can be derived as [26] 

eiP, X) = [^]F{X) + (^)(1 - Pp-) (1) 

where kp — (l.STr^p)^ which is equal to Fermi momentum in case of SNM, the 
kinetic energy per nucleon e'^^" = [^^]F{X) with F{X) = [(i±^)!^!+(lz2Q!f!] 
and = J^oQ + X'^J^Qi where J^,qo and J^qi represent the volume integrals of 
the isoscalar and the isovector parts of the M3Y interaction supplemented by 
the zero-range potentials having the form 



Jvoo ^ Jvooie'n - J J jC{s,e)d's 

A-jr A-jT 

= 7999— - 2134— J + Joo(l - ae'^'") where Jqo = -276 MeV fm^ (2) 



Jvoi ^ Jvoiie'n - 1 1 jC'{s,e)dh 
4:71 Att 

= -4886— + 1176 — - + Joi(l - ae''''') where Jqi = 228 MeV fm^ (3) 
4'^ 2.5'^ 

The isoscalar t^g^^ and the isovector tgi'^^ components of M3Y interaction po- 
tential [19,26] supplemented by zero range potentials are given by t^^'^is, e) — 
79992iPl_M_2l342^PgM_276(l-ae)<5(5) andt^^y^^^^) ^ _48862^4j^ + 
^ 228(1 - ae)S{s) respectively, where the energy dependence 
parameter a=0.005/MeV. The DDM3Y effective NN interaction is given by 
VQi{s,p,e) = t^j^^^ {s,e)g{p) where the density dependence g{p) = C(l — /3p") 
and the constants C and f3 of density dependence have been obtained from 
the saturation condition |^ = at p = po and e = eo where po and eo are 
the saturation density and the saturation energy per nucleon respectively. 
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The equilibrium density of the cold SNM is determined from the saturation 
condition. Then Eq.(l) along with the equation for the saturation condition 
mentioned above can be solved simultaneously for fixed values of the satu- 
ration energy per nucleon eo and the saturation density po of the cold SNM 
to obtain the values of /3 and C. The constants of density dependence /3 and 

C, thus obtained, are given by 3 — , / ''l^" -iom where p = tt^tt4, 

q = with J° = J^,oo(eo^") which means J^oo at 6^=*" = eg^", the kinetic 
energy part of the saturation energy per nucleon of SNM, kpo = [l.Svr^po]^^^ 
and C = ^° gn^k^ (i-^p") respectively. It is quite obvious that 

5mJgooPo[l-(n+l)/3pg- ' fg^,^ ] 

the constants of density dependence C and f3 obtained by this method depend 
on the saturation energy per nucleon cq, the saturation density po, the index 
n of the density dependent part and on the strengths of the M3Y interaction 
through the volume integral J°oo- 

The incompressibility or the compression modulus of the SNM, which is a 
measure of the curvature of an EoS at saturation density and defined as 
kF§j^ \kF=kFQ^ measures the stiffness of an EoS can be obtained theoretically: 



dl4 



\kF=kFn 



dp^ 



P=Pa 



Sh'k% 

5m 



9J^,,Cn{n + l)Pp^ 



,n+l 



-9aJooC[l - (n + 1)(3pq\— ^ H S- 



5m 



10m 



(4) 



The EoS of isospin asymmetric nuclear matter, given by Eq.(l) can be gener- 
ally expanded as 

e(p, X) = e(p, 0) + Esym{p)X^ + (5) 



and Esyraip) = \ ^ ^qx^^ \x=o is the nuclear symmetry energy. The absence 
of odd-order terms in X in Eq.(5) is due to the exchange symmetry between 
protons and neutrons in nuclear matter when one neglects the Coulomb inter- 
action and assumes the charge symmetry of nuclear forces. The higher-order 
terms in X are negligible and as a good approximation, the density-dependent 
symmetry energy Esymip) can be extracted from [27] 

E,ym{p)^e(p,l)-e(p,0) (6) 



which can be obtained using Eq.(l) and represents a penalty levied on the 
system as it departs from the symmetric limit of equal number of protons and 
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neutrons and can be defined as the energy required per nucleon to change the 
SNM to pure neutron matter (PNM). 



3 Isospin dependent properties of asymmetric nuclear matter 



Around the nuclear matter saturation density po the nuclear symmetry energy 
Esyra{p) can be expanded to second order in density as 



where L and Ksym are the slope and curvature parameters of nuclear symmetry 
energy at po a-^d hence 

^ — 9p — ^"^'^ ^ ' 



K - 9o2^^^£y»n(p) I /Q^ 

^^sym — ^Fo \p=PO \'^J 



The L and Ksym characterize the density dependence of the nuclear symme- 
try energy around normal nuclear matter density and thus carry important 
information on the properties of nuclear symmetry energy at both high and 
low densities. In particular, the slope parameter L has been found to corre- 
late linearly with the neutron-skin thickness of heavy nuclei and thus can in 
principle be determined from the measured thickness of neutron skin of such 
nuclei [5,9,28,29,30,31,32,33]. Recently, this has been possible [24] although 
there are large uncertainties in the experimental measurements. 

The Eq.(6) can be differentiated twice with respect to nucleonic density p 
using Eq.(l) to yield 



op o P Po ^ 

-^4^E°(^)2/3[l-/?p"]F(l) 
5 Po 

5 Po 

-:^(22/3 _ V)aJ^^CEl{^fl\l - (n + l)/3p1 (10) 
10 Po 
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9p2 15' Po 2 

2q;JoiCE^^ p_ 2/3 



5 p Po' 
io p Po 

• it) p 



(22/3_l)^Ml:^(^)2/3[i_^^n] 

15 p Po 



-^(2^^ - l)aJ^C^{^fl'[l - (n + l)/3p"] 
5 p Po 

+ 3 (2^/3 - l)aJooCE%{-^flMn + l)/3p"-i (11) 
iU Po 

where the Fermi energy E^p = for the SNM at ground state. The above 
Eqs.(10,ll) at p=po are used to evaluate the values of L and Kgym using 
Eqs.(8,9). At the nuclear matter saturation density and around X = 0, the 
isobaric incomprcssibility of asymmetric nuclear matter can also be expressed 
to second order in X as [34,35] KiX^ ^ fCoo + KasyX'^ where Koo is the 
incomprcssibility of symmetric nuclear matter at the nuclear matter saturation 
density, and the isospin dependent part [36] 

^asy ~ ^sym ~ 6L (12) 



characterizes the density dependence of the nuclear symmetry energy. Infor- 
mation on Kasy can in principle be extracted experimentally by measuring the 
GMR in neutron-rich nuclei. 



4 Results of present calculations and experimental status 



The calculations are performed using the values of the saturation density 
Po=0.1533 fm~3 [37] and the saturation energy per nucleon eo=-15.26 MeV 
[38] for the SNM obtained from the co-efficient of the volume term of Bethe- 
Weizsacker mass formula [39,40] which is evaluated by fitting the recent ex- 
perimental and estimated atomic mass excesses from Audi-Wapstra-Thibault 
atomic mass table [41] by minimizing the mean square deviation incorpo- 
rating correction for the electronic binding energy [42]. In a similar recent 
work, including surface symmetry energy term, Wigner term, shell correction 
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and proton form factor correction to Coulomb energy also, turns out to 
be 15.4496 MeV [43] (a^ =14.8497 MeV when A° and A^/^ terms are also 
included). Using the usual values of q;=0.005 MeV~^ for the parameter of en- 
ergy dependence of the zero range potential and n—2/3, the values obtained for 
the constants of density dependence C and P and the SNM incompressibihty 
Koo are 2.2497, 1.5934 fm^ and 274.7 McV respectively. The saturation en- 
ergy per nucleon is the volume energy coefficient and the value of -15.26±0.52 
MeV covers, more or less, the entire range of values obtained for for which 
now the values of C=2.2497±0.0420, /9=1.5934±0.0085 fm^, the SNM incom- 
pressibihty K^—274:.7±7.4: MeV, nuclear symmetry energy at nuclear density 
^sj/m(po)=30.71 ±0.26 MeV, the slope L=45.11±0.02 MeV and the curvature 
-f^sym=-183.7±3.6 MeV parameters of the nuclear symmetry energy and the 
isospin dependent part ii'asj^=-454.4±3.5 MeV of the isobaric incompressibil- 
ity. These values are compared with those obtained from few other theoretical 
models in Table 1. In Fig.l the NSE is plotted as a function of p/po for the 
present calculation using DDM3Y interaction, and compared with those for 
Akmal-Pandharipande-Ravenhall [44] and MDI interaction [45]. 

The theoretical estimate of K^^, from the refractive a-nucleus scattering is 
about 240-270 MeV [46,47]. The recent experimental determination of Koo 
based upon the production of hard photons in heavy ion collisions which led 
to the experimental estimate of = 290 ± 50 McV [48] . However, the ex- 
perimental values of K^o extracted from the isoscalar giant dipole resonance 
(ISGDR) are claimed to be smaller [49]. Considering the status of experimen- 
tal determination of the SNM incompressibihty from data on the compression 
modes ISGMR and ISGDR of nuclei it can be inferred [50] that due to viola- 
tions of self consistency in HF-RPA calculations of the strength functions of 
giant resonances result in shifts in the calculated values of the centroid energies 
which may be larger in magnitude than the current experimental uncertain- 
ties. In fact, the prediction of lying in the range of 210-220 MeV were due 
to the use of a not fully self-consistent Skyrme calculations [50]. Correcting 
for this drawback, Skyrme parmetrizations of SLy4 type predict Koo values 
in the range of 230-240 MeV [50]. Moreover, it is possible to build bona fide 
Skyrme forces so that the SNM incompressibihty is close to the relativistic 
value, namely 250-270 MeV. Therefore, from the ISGMR experimental data 
the conclusion can be drawn that ^ 240 ± 20 MeV. The ISGDR data tend 
to point to lower values [49,51,52] for Koo- However, there is consensus that the 
extraction of Koo is in this case more problematic for various reasons. In par- 
ticular, the maximum cross-section for ISGDR decreases very strongly at high 
excitation energy and may drop below the current experimental sensitivity 
for excitation energies [50] above 30 and 26 MeV for ^^^Sn and ^°^Pb, respec- 
tively. The present non-relativistic mean field model estimate for the nuclear 
incompressibihty Koo for SNM using DDM3Y interaction is rather close to the 
theoretical estimates obtained using relativistic mean field models and close 
to the upper limit of the recent values [53] extracted from experiments. The 
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Fig. 1. The NSE (nuclear symmetry energy Egym) is plotted as a function of p/po for 
the present calculation using DDM3Y interaction and its comparison, with those for 
Akmal-Pandharipande-Ravenhall (APR) [44] and MDI interaction for the variable 
x=0.0, 0.5 defined in Ref. [45]. 

present consensus is that the acceptable value [54,55] for the incompressibility 
of symmetric nuclear matter lies in the range of 250-270 MeV. The present 
value of 274.7±7.4 MeV for the incompressibility K^o of SNM obtained using 
DDM3Y interaction is, therefore, an excellent theoretical result. 

The constant of density dependence /5=1.5934±0.0085 fm^, which has the di- 
mension of cross section for n=2/3, can be interpreted as the isospin averaged 
effective nucleon-nucleon interaction cross section in ground state symmetric 
nuclear medium. For a nucleon in ground state nuclear matter kp ^ 1.3 fm~^ 
and go ~ hkpc ~ 260 MeV and the present result for the 'in medium' effec- 
tive cross section is reasonably close to the value obtained from a rigorous 
Dirac-Brueckner-Hartree-Fock [56] calculations corresponding to such kp and 
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qo values which is ^ 12 mb. Using the vahic of the constant of density de- 
pendence /?=1.5934±0.0085 fm^ corresponding to the standard value of the 
parameter n=2/3 along with the nucleonic density of 0.1533 fm~^, the value 
obtained for the nuclear mean free path A is about 4 fm which is in excellent 
agreement [57] with that obtained using another method. 

The volume symmetry energy coefficient extracted from the masses of fi- 
nite nuclei provides a constraint on the nuclear symmetry energy at nuclear 
density EgymiPo)- The value of = 30.048 ± 0.004 MeV recently extracted 
[58] from the measured atomic mass excesses of 2228 nuclei is reasonably 
close to the theoretical estimate of the value of NSE at the saturation den- 
sity -Esj/m(po)=30.71±0.26 MeV obtained from the present calculations using 
DDM3Y interaction. In ref. [59] it is between 29.10 MeV to 32.67 MeV and 
that obtained by the hquid droplet model calculation of ref. [3] is 27.3 MeV 
whereas in ref. [60] it is 28.0 MeV. It should be mentioned that the value of 
the volume symmetry parameter Sy in some advanced mass description [61] is 
close to the present value which with their —KvoI-KoI = equals 29.3 MeV. 
The value of NSE at nuclear saturation density ~ 30 MeV, therefore, seems 
well established empirically. Theoretically different parameterizations of the 
relativistic mean-field (RMF) models, which fit observables for isospin sym- 
metric nuclei well, lead to a relatively wide range of predictions 24-40 MeV for 
Esymipo)- The present result of 30.71±0.26 MeV of the mean field calculation 
is close to the results of the calculation using Skyrmc interaction SkMP (29.9 
MeV) [62], Avl8+(5t;+UIX* variational calculation (30.1 MeV) [44] and field 
theoretical calculation DD-F (31.6 MeV) [27]. 



Information on K^sy can in principle be extracted experimentally by measuring 
the GMR in neutron-rich nuclei. Earlier attempts based on this method have 
given, however, widely different values. A value of Kasy = — 320±180 MeV [69] 
with a large uncertainty was obtained from a systematic study of the GMR in 
the isotopic chains of Sn and Sm. In this analysis, the value of Koo was found to 
be 300±25 MeV. The recently obtained values from the isoscaling analyses of 
isotope ratios in intermediate energy heavy-ion collisions [22] gives L Ri65 MeV 
and Kasy ~ -453 MeV. The extracted values of L = 88±25 MeV and K^sy = 
-500±50 MeV from the isospin diffusion data [70] and Kasy = -550±100 
MeV obtained from recently measured isotopic dependence of the GMR in 
even-A Sn isotopes [23] are consistent with the value of K^sy obtained from 
the present calculations. The incompressibility of symmetric nuclear matter 
at its saturation density has been determined to be 240±20 MeV [71] from 
analyses of the GMR. A very recent extraction from neutron skin thickness of 
nuclei provides the values of L=75±25 MeV (revised to narrower window of 
L-45-75 MeV recently [72]) and Kasy=-500tlll MeV [24] and it is needless 
to stress that the theoretical values of our calculations for L=45.11±0.02 MeV 
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Table 1 

Results of the present calculations (DDM3Y) of incompressibility of isospin symmet- 
ric nuclear matter K^, nuclear symmetry energy at saturation density Esym{po), 
the slope L and the curvature Kaym parameters of the nuclear symmetry energy 
and the isospin dependent part Kasy of the isobaric incompressibility (all in MeV) 
are compared with those obtained with other models. 



Model 




EsymiPo) 


L 




Kasy 


Ref. 


Expt. 





31.6 


75 ±25 







[24] 


This work 


274.7 ± 7.4 


30.71 ± 0.26 


45.11 ±0.02 


-183.7 ±3.6 


-454.4 ± 3.5 




NLl 


212 


43.5 


140 


143 


-697 


[63] 


NL2 


401 


44.0 


130 


20 


-750 


[63] 


NL3 


271 


37.3 


118 


100 


-608 


[64] 


NL-SH 


356 


36.1 


114 


80 


-604 


[65] 


FSUGold 


229 


32.5 


60 


-52 


-412 


[5] 


DD-MEl 


245 


33.1 


55 


-101 


-431 


[66] 


DD-ME2 


251 


32.3 


51 


-87 


-393 


[67] 


DD 


241 


31.7 


56 


-95 


-431 


[68] 


DD-F 


223 


31.6 


56 


-140 


-476 


[27] 



and Kasy— —4:54.4 ± 3.5 MeV are in excellent agreement. 



5 Summciry and conclusion 

We show that the theoretical description of nuclear matter based on mean 
field calculation using DDM3Y effective NN interaction yields a value of nu- 
clear incompressibility which is highly in agreement with that extracted from 
experiment and gives a value of symmetry energy that is consistent with the 
empirical value extracted by fitting the droplet model to the measured atomic 
mass excesses and with other modern theoretical descriptions of nuclear mat- 
ter. The present NSE is 'soft' because it increases initially with nucleonic 
density up to about two times the normal nuclear density and then decreases 
monotonically at higher densities and is consistent with the recent evidence 
for a soft NSE at suprasaturation densities [45]. This interaction also pro- 
vides good descriptions for the elastic and inelastic scattering [73] and various 
kinds of radioactivities [74] . The slope L and the isospin dependent part Kasy 
of the isobaric incompressibility are consistent with the constraints recently 
extracted from analyses of experimental data. Of all other models, DD-ME2 
provides comparatively better estimates and of all RMF models provides best 
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mass predictions with r.m.s. deviation of 0.9 MeV. Interestingly, our calcula- 
tions provide so far the best theoretical estimates for the isospin dependent 
properties of asymmetric nuclear matter. 
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